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ABSTRACT 



Various computational methods and operational computer codes used to predict the 
aerodynamic coefficients and separation trajectories of aircraft stores are examined. The 
semi-empirical aeroprediction code Missile DATCOM is used to obtain the coefficients 
of a modeled store. These coefficients, together with the modeled ejection forces, are 
used in free-stream state-space equations of motion to predict the store trajectory. The 
results are compared with the Nielson Engineering and Research (NEAR) store separation 
code which provides accurate trajectory profiles, for speeds below the critical speed, by 
use of a vortex-lattice and panel method. Modification of the Missile DATCOM 
aerodynamic coefficients provides single-point state-space prediction of the store pitch 
trajectory within 30% of the NEAR code results. Store trajectories were restricted to the 
first 0.2 seconds of free flight. 
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STORE SEPARATION INTRODUCTION 



I . 

A . BACKGROUND 

The prediction of the trajectory of a store ejected from 
an aircraft is of major concern to defense aviation. The 
increasing requirement for aircraft to fulfill a multiple role 
in the hostile environment demands that it also carry an ever- 
increasing variety of stores . These stores range from 
missiles and bombs stockpiled for many years to new, 
aerodynamically complex stores with lifting bodies. 

The aim of this work is to review current methods of 
predicting store releases, compare results of two different 
methods, and to present an easy-to-use methodology for 
investigating the subsonic release of a unsophisticated store. 

A non-axisymmetric pod ejected from the F/A-18 outboard 
pylon is presented as an example. This store was chosen 
because it is of practical value due to the extensive use of 
these modified pods at the Pacific Missile Test Center, at 
Point Mugu, California. Non-axisymmetrical stores are very 
difficult to model and present a special problem in the store 
separation field. This pod is, therefore, modeled as 
axisymmetric and non-axisymmetric, and the results are 
compared. In addition, the store lacks aerodynamic control 
surfaces and is highly unstable. Application of the 
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methodology presented here should eventually provide the user 
with accurate results for a minimal time and cost expenditure. 

B . MATHEMATICAL MODELING 

There are many ways to mathematically model store 
separation. Economics and accuracy are always the main 
considerations when exploring the codes and methods to be 
used. For simple cases, intermediate approaches such as panel 
methods and solutions of the Euler equations provide 
sufficient accuracy. Calculations involving flow separation 
and shock wave interference become too complex for 
intermediate methods and acceptable results can only be 
obtained using modern computational fluid dynamic (CFD) 
techniques involving solutions of the Navier-Stokes equations 
[Ref. l:p. 7-2] . For speeds below the subsonic Mach critical 
speed, component buildup methods (empirical and semi- 
empirical) , and panel methods provide the necessary accuracy 
and are much easier to use. The codes chosen for this study 
are explained in detail in their respective sections. 

The F/A-18 model was obtained from Mr. L. L. Gleason of 
the Ordnance Systems Department, Naval Weapons Center, at 
China Lake, California. The F/A-18 was chosen due to its wide 
variety of store carriage and also because of its solid future 
with Naval aviation. 

The store model is based upon both the Missile Datcom 
method [Ref 2] and also by the Nielson program [Ref 3] . These 
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will be discussed in Sections III and IV. Appendix A contains 
the physical description of the pod and its inertial 
characteristics . 

C . METHODOLOGY 

The initial investigation (Section II) consisted of 
looking at the ejector forces, combined with the store's 
inertial characteristics, to give a sense of the magnitude of 
initial velocities and moments. These forces and moments, 
while easy to predict, were insufficient to provide real 
evidence of a safe release or of a release problem. A large 
pitch-down moment or a large vertical velocity may be off-set 
by unforeseen aerodynamic forces . The velocities and moments 
calculated will be used in Section III and compared with 
results of Section IV for trajectory prediction. 

Section III introduces a method of predicting the store's 
longitudinal and lateral aerodynamic coefficients using a 
missile code (Missile DATCOM) developed by McDonnell Douglas 
Missile Systems Company for the Flight Dynamics Laboratory 
(FDL) , Wright-Patterson Air Force Base (WFAFB) , Ohio [Ref. 2] . 
These aerodynamic coefficients are then used in longitudinal 
and lateral equations of motion to predict the free-stream 
trajectory of the store. A brief explanation of the 
derivation of these equations of motion is provided to the 
reader. An accurate representation of the aerodynamic 
characteristics is essential for the full prediction of the 
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shape's trajectory. Missile DATCOM is based on the body 
buildup method and includes a number of prediction methods for 
each component of the configuration. Other codes, such as 
MISL3, are also available, each with its own relative 
strengths . 

Section IV describes a computer prediction method 
developed by Nielson Engineering and Research, Inc., (NEAR), 
also under contract to the United States Air Force. This code 
provides a six degree-of-f reedom (6D0F) simulation which takes 
into account the aircraft flowfield using vortex-lattice and 
panel methods. The NEAR simulation provides a high fidelity 
representation for the subsonic case, but also requires the 
most effort and computer capability. An attempt has been made 
to minimize the complexity of the input to the program. Once 
a good data base has been developed the NEAR code should prove 
easy to use. It has been used extensively throughout the 
defense industry, and has undergone many modifications to 
incorporate improvements and options . 

Section V compares the trajectory results of the 
linearized aerodynamic simulation method and the NEAR code. 

A modification to the linear coefficients used in the linear 
aerodynamic simulation is also discussed. 
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II. EJECTION FORCES AND MOMENTS 



Although quite uncomplicated in nature, calculating the 
ejection forces and moments provides the user with a feel for 
the magnitude of the initial movement of the store. This 
procedure should be done prior to using any other method as a 
preliminary step in order to prepare for the simulation. 

A. EJECTOR CHARACTERISTICS 

The non-aerodynamic forces involved in the ejection of a 
store depend upon the ejector cartridge, the ejector rack, the 
weight of the store, and the rigidity of the wing. Each of 
these parameters determine the resultant moments and forces 
provided to the store. The store and the wing (or fuselage) 
were considered as rigid, thus neglecting any aeroelastic or 
structural bending effects. None of the methods discussed 
here address aeroelastic bending due to the added complexity 
of the problem. 

The Douglas production BRU-32/A Bomb Ejector Rack combines 
two sets of hooks, one set at 14-inch spacing and one set at 
30-inch spacing, with an ejection system designed for carriage 
of stores with suspension lugs per MIL-A-8591. Two 
electrically initiated CCU-45B cartridges are used for store 
release and ejection. The self-retracting ejector pistons, 
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spaced symmetrically 20 inches apart at each end of the rack, 
have a piston stroke of six inches. These pistons are spring 
loaded against the store during loading to prevent impact of 
pistons during firing. The orifice sizes can be varied, by 
replacement, to provide force and pitch control for store 
separation. [Ref. 4] 

The best source of ejector data is from the Aircraft 
Ordnance Procedures (AOP) contained in the Aircraft Stores 
Interface Manual, (Reference 4) . This manual contains a 
complete description of the rack under consideration. Here, 
the example separation uses two 14-inch spaced hooks and 

0. 118-Diameter orifices. The force diagram corresponding to 
this orifice, modeled from the AOP, is shown in Figure 1. 
Sufficient data points were read from the AOP graph in order 
to curve-fit the points using Computer Associates' 
CricketGraph® software [Ref. 5]. This provided two things; 

1. ) integration of the polynomial equation gave the total 
impulse value, and 2.) Nielson software in Section IV uses a 
fifth-order polynomial for calculations. Note that the AOP 
shows the total force provided by the ejector, while the NEAR 
simulation requires a f orce-per-e jector-f oot polynomial. 
Given the shape of the force curve, a triangular 
representation of the curve gave a value of 348 lb f -sec versus 
the integrated value of 344.9 lb f -sec, which corresponds to an 
error of only 1%. Therefore elementary estimates of the total 
impulse are practical. 



6 



FORCE (LBF) 




0.000 0.010 0.020 0.030 0.040 0.050 

TIME (SECS) 



Figure 1 Ejector Force 
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FORCES 



B. 

1 . Theory 

Using the conservation of linear momentum theory 
presented in Appendix B, the pod was modeled as a simple beam. 
The weight was concentrated at the center of gravity and the 
resultant ejection forces were placed at the ejector feet 
locations. The ejection force time histories are treated as 
an equivalent force impulse which results in a change in 
linear momentum as shown by the following equation: 

JZF, dt = G 2 - G x = (M*V t ) 2 ~ (M*V Z ) x 

The integral for representing the linear impulse is 
based upon the fifth-order polynomial force time history shown 
here . 

y= 280. 2-608022. 8x + 180688545 . 8x 2 - 10638451765 . 5x 3 + 
234526098234 . 3x 4 - 1783873760363 . 4x 5 

2 . Example 

For our example, the BRU-32A bomb rack used two CCU- 
45B cartridge-activated devices (CADS) . The peak force was 
14,500 lb f and the pulse duration was 50 msec. The total 
linear impulse calculated was 345 lb f -sec . Using a weight of 
371 pounds for the pod, the end-of-stroke velocity was 29.9 
feet per sec (fps) . This value was used in Section III as the 
initial velocity for the trajectory simulation. 
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The graphs contained in the AOP shows end-of-stroke 
velocity for a 371 lb store to be approximately 24.5 fps . 
This graph is derived from empirical data and should provide 
the most accurate representation of the actual release 
velocity. The values of 24.5 and 29.9 fps are compared with 
the end-of-stroke store velocity of 30.9 fps predicted from 
the NEAR code in Section IV. 

C . MOMENTS 

1 . Theory 

The pod was modeled as a simple beam. The inertial 
characteristics of the store are listed in Appendix A. The 
two ejector feet provide the beam with a moment. The 
equations used are elementary, and with certain assumptions 
the angular velocity can be calculated. The equations are 
based on conservation of angular momentum. The equations and 
assumptions are presented in Appendix B for completeness. The 
resulting equation is shown here. 

jlMy dt = Iyy * (By 

2 . Example 

The center of gravity was offset aft of the center of 
the ejector feet by 4.05 inches. The impulse thus provided 
the pod with 1.12 radians per second (rps) or 64.1 deg/sec 
initial angular velocity. The question now is, are these 
reasonable values? At the end of 0.429 seconds, the shape 
will be one body-length below the aircraft. At the same time, 
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the shape will have rotated 27.5 degrees nosedown. Therefore, 
without considering the aerodynamics of the vehicle, the pod 
seems to have cleared the wing. A good rule of thumb is 2-3 
body lengths clearance, although exceptions do occur. 



D . PRELIMINARY OBSERVATION 

This initial calculation provides an approximate estimate 
of the forces and moments involved. The values obtained are 
of sufficient accuracy for calculations such as end-of-stroke 
loading, etc. Excessively large separation velocities are 
usually an indication of a miscalculation. Translational 
velocities should range from 10-30 fps and pitching velocities 
from 0-6 rps . 
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III. AERODYNAMIC CALCULATION 



A , THEORY 

Many physical systems can be modeled by second-order 
differential equations. The mathematical treatment of fixed- 
wing flight vehicle motions was first developed by G.H. Bryan. 
He laid the mathematical foundation for airplane dynamic 
stability analysis, developed the concept of the aerodynamic 
stability derivative, and recognized that the equations of 
motion could be separated into a symmetric longitudinal motion 
and an unsymmetric lateral motion. Experimental studies were 
initiated by L. Bairstow and B.M. Jones of the National 
Physical Laboratory in England, and by Jerome Hunsaker of the 
Massachusetts Institute of Technology to determine estimates 
of the aerodynamic stability derivatives in Bryan's theory. 
In addition to determining stability derivatives from wind- 
tunnel tests of scale models, Bairstow and Jones 

nondimensionalized the equations of motion and showed that, 
with certain assumptions, there were two independent 

solutions, i.e., one longitudinal and one lateral. These two 
solutions provide the free-stream trajectories we seek. [Ref. 
6 : p . 113] 

The dynamic stability characteristics of a pod can be 
represented by six equations of motion, three for the forces 
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involved X, Y, and Z and three for the moments L, M and N. 
The force equations relate the forces acting on the body to 
the corresponding linear accelerations and the moment 
equations relate the moments to the corresponding angular 
accelerations. It is usually possible to consider the 

longitudinal motions completely separately from the lateral- 
directional motion, by neglecting the various coupling terms. 
[Ref. 7 : p . 14] 

As a caveat to the use of this method, these equations are 
the linearized version and are only valid up to approximately 
10 degrees angle-of-attack . Treatment of the non-linear 
aerodynamic coefficients, while not extremely difficult, does 
require knowledge of the behavior of the coefficients. Since 
the goal of this report is to predict the behavior of new 
shape configurations, this knowledge is not presumed. 
Therefore, the results obtained in this fashion are only valid 
for estimating the motion in the first fractions of a second. 
These results are compared with those obtained in the more 
rigorous method of Section IV. 

1 . Longitudinal Equations 

The rigid-body, longitudinal equations of motion can 
be developed from Newton's second law: 

^Pitching moments=ZM cg =I yy *0 

The assumption that body motion consists of small deviations 
from its equilibrium flight condition allows us to use 
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perturbation theory to examine the aerodynamic force 
derivatives in terms of angle-of-attack (AOA) , vertical 
velocity, pitch angle, and pitch rate, by means of a Taylor 
series expansion. The X-force, Z-force, and pitching moment 
equations comprise the longitudinal equations. To separate 
these equations from the lateral equations, they must not be 
coupled. This is a reasonable assumption given the nominal 
geometric shape of missiles or pods [Ref. 6] . The body is 
constrained to move in a vertical plane and is free to pitch 
about its center of gravity. For a comprehensive derivation 
of the equations of motion, see Reference 6. The resulting 
equations of motion are shown in Appendix C. 

2 . Lateral Equations 

The Y-force, rolling, and yawing moment equations 
comprise the lateral equations of motion. Once again, these 
equations are uncoupled from the longitudinal equations by the 
assumption of small cross-component moments of inertia. That 
is, I zy , I xt , I yx are small in comparison to the principal axis 
moments. These equations are also valid assuming only small 
variations in displacement and velocity. The lateral 
equations are derived from the following Newton's laws: 
ZRolling moment s= I „*<!> 

ZYawing moments=I Il *\|/ 

The third equation is derived from the Taylor series expansion 
of the side force derivative. The lateral motion of the pod 
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disturbed from its equilibrium state is a complicated 
combination of rolling, yawing, and sideslipping motions. 
Assuming the cross products of inertia are ignored, some of 
the coupling terms can be simplified. 

3 . State Variable Representation 

The linearized longitudinal and lateral equations 
developed above are simple, ordinary linear differential 
equations with constant coefficients. The coefficients in the 
differential equations are made up of the aerodynamic 
stability derivatives, mass, and inertia characteristics of 
the pod. These equations can be written as a set of first- 
order differential equations in the state-space (state 
variable) form: 



{k} = [A] * { x) + [B] * {T| } 

where {x} is the state vector, {T|} is the control vector and 
the matrices [A] and [B] contain the pod' s dimensional 
stability derivatives. In our case the pod does not have any 
active control surfaces. For missile launches, however, the 
control vector and the B-matrix would be used to represent the 
control surfaces and control input. The simplified state- 
space representations of the longitudinal and lateral 
equations of motion are shown in Appendix C. The required 
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coefficients listed in Appendix C are derived using the 
procedures in Section III.B. 

The state-space representation of the equations of 
motion can be solved simultaneously using matrix software such 
as MATLAB® [Ref. 8] or Control-C® [Ref. 9] . Given the 
conditions of the ejection, namely, initial translational 
velocities, initial angular velocities, and attitude angles, 
solutions of the state-space equations can be used to predict 
the trajectory of the pod. The initial angles and velocities 
are the initial conditions imposed upon the state-space 
equations. The Control-C® commands for this procedure are 
contained in Appendix D. 

B. MISDATCOM 

MISDATCOM was developed by McDonnell Douglas Astronautics 
Company, St Louis, Missouri, for the Flight Dynamics 
Laboratory of the Air Force Wright Aeronautical Laboratories, 
Wright -Patterson AFB . The program was completed in December 
1985. The Missile Datcom was created to provide an 
aerodynamic design tool which has the predictive accuracy 
suitable for preliminary design, and the capability for the 
user to easily substitute methods to fit specific 
applications. [Ref. 2:p. 1] 

1 . Theory 

There are many different methods to predict the 
aerodynamic static, dynamic, and control characteristics of 
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missiles. Component build-up was chosen as the most suitable 
for this program. Although panel methods are better suited 
for arbitrary configurations, component build-up was chosen 
due to the accuracy provided for conventional configurations 
and for the ability to do parametric studies easier. 

Basically, component build-up consists of using 
various methods to compute the characteristics (skin friction, 
force and moment coefficients, panel loading, . . .) of the 
individual configuration components. The various methods are 
chosen for their applicability to the configuration or flight 
condition. Then the components are combined. Previous 
methods of combining the components (fins, body, engine inlet) 
consisted of adding the individual coefficients and then 
multiplying the sum by some interference factor obtained using 
slender body theory. The approach taken with MISDATCOM was to 
use the "equivalent angle of attack method" developed by 
Nielson Engineering and Research, Inc. (NEAR). This method 
assumes that the panel loading for a given panel angle-of- 
attack is unique. With this method the panel angle of attack 
is computed including the effect of panel roll orientation 
with respect to the free stream velocity vector, panel 
proximity to the fuselage or to other panels, and external 
vortex flow field effects. Then the isolated panel 
characteristics are interpolated at the panel equivalent angle 
of attack to yield the panel load when mounted on a body in 
combination with other surfaces. [Ref. 2] 
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2 . Procedure 



The procedure for using MISDATCOM is straightforward. 
The pod shape is modeled using simple geometric shapes. The 
previous release of MISDATCOM software could only handle 
axisynunetric shapes or forms with a vertical plane of 
symmetry. Unfortunately, the pod under consideration is non- 
axisymmetric and therefore the results are not entirely valid. 
The latest version, however, can model some non-axisymmetry 
through the use of a "protuberance" option. This is the April 
1991 release and is now available. 

For this investigation, a simulation run was made 
using both versions and a comparison of the aerodynamic 
coefficients was made. The pod radius was also varied from 
the minimum pod to the maximum, and the resulting coefficients 
compared. The MISDATCOM code coefficients were somewhat 
insensitive to radius changes of this magnitude, therefore, 
the minimum-radius, axisymmetric case was retained for 
comparison with the minimum-radius non-axi symmetric case. 
Figure 2 shows the pod semi-profile of the MISDATCOM procedure 
alongside the pod semi-profile of the NEAR simulation. Figure 
3 contains a sample input to the program and Figure 4 lists a 
partial sample output . The actual input and output of the 
program are contained in Appendix E. These aerodynamic 
coefficients are next used in the longitudinal and lateral 
state-space equations of motion to predict the store 
trajectory . 
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Figure 4 MISDATCOM Output 
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C . AERODYNAMIC CALCULATION 



1 . Longitudinal 

The aerodynamic coefficients derived from the 
MISDATCOM were entered into the plant matrix [A] of the 
longitudinal equations of motion. The initial conditions were 
applied and time history of the pitch angle, pitch angular 
velocity, and vertical velocity was found. Figures 5-8 show 
that the pod assumes a large nose-down, highly divergent 
pitching motion. This is understandable since there are no 
control surfaces to make the pod stable. Examination of the 
roots of the plant matrix [A] indicates an unstable flight 
vehicle. The pitch angle and pitch rate shown are in the 
store body coordinate frame. 

These results are valid only for the range of linear 
values of the aerodynamic coefficients . This is approximately 
up to 10-12 degrees angle-of-attack . Therefore, an estimate 
of the non-linear behavior of the pod after an AOA of 10 
degrees is reached is necessary. Due to the large pitch 
angular velocity, it is obvious that this limit is reached 
after only 0.2 seconds. Section V discusses the acceptable 
time range where these results are valid. 

To estimate the non-linear behavior, the dominating 
terms in the [A] -matrix must be determined. This provides an 
insight into the possible range of values to substitute into 
the [A] -matrix. A discussion of a possible approach to this 
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Figure 5 Pitch Angle Response 




Figure 6 Pitch Rate Response 
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Figure 7 Vertical Separation Distance 




Figure 8 



Vertical Velocity 
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problem is presented in Section V. To completely investigate 
this area is beyond the scope of the method outlined here. 
Comparison of the "linear" trajectory obtained here will be 
made with the trajectory obtained in Section IV, which does 
take into account the non-linearity of the coefficients. 

2 . Lateral 

Due to the nearly vertical forces and moments provided 
by the ejection rack during straight and level flight, the 
free-stream lateral equations of motion will not provide us 
with any insight into the safe jettison of the shape. 
However, if the pod were experiencing sideslip, then the 
sideslip could be entered into the state-space equations as an 
initial condition. For this pod a sideslip of -0.5 degrees 
was assumed. The resulting lateral motion is not shown 
because the main emphasis is on the longitudinal motion. The 
pod is unstable laterally, also, but the initial movement is 
small due to the relatively small initial conditions. In 
addition to the yaw angle and rate, if the ejection rack 
provided an initial roll rate, such as when the pod is loaded 
off-center of the bomb rack, then that influence could also be 
included . 

D. PARAMETER VARIATION 

The parameter variation due to non-axisymmetry is 
difficult to handle. Comparison with shapes with known 
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(experimentally obtained) coefficients might provide some 
degree of accuracy in predicting more accurate results . 

E. PRELIMINARY OBSERVATION 

The results of the method in Section II were a vertical 
velocity of 29.9 feet/sec and an angular velocity of 64.1 
deg/sec. The aerodynamic method of Section III begins with a 
vertical velocity of 29.9 fps, quickly diverging to an 
extreme value. The pitch rate, q, also diverges quickly. 
Because the pod is very unstable any linear approximation of 
its behavior will have strict limits. The coefficients used 
in this calculation, from the MISDATCOM code, are not valid 
beyond a fraction of the trajectory. The valid time range of 
prediction is presented in Section V. 
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IV. 



"NEAR" PROGRAM 



A . BACKGROUND 

A computer prediction method was developed by Nielsen 
Engineering and Research, Inc., (NEAR) under contract to the 
United States Air Force. The work was performed during the 
period 1968 to 1972. The final result is a method for 
predicting the six degree-of-f reedom store separation 
trajectory at speeds up to the subsonic critical Mach number. 
After delivery of the program to the Government, many new 
capabilities were added. The code used for this paper was 
obtained from the Naval Weapons Center, China Lake, Ca . The 
program has been widely accepted by industry and government . 
The code encompasses 9, 900 lines of code and thus requires a 
device with sufficient computer memory for operation. (This 
also depends upon the program application.) 

The aircraft fuselage, separated store body, and adjacent 
stores are modeled using point sources and sinks. Angle of 
attack effects are included using a cross-flow model. The 
aircraft wing and wing pylons are modeled using planar vortex 
lattice models which include dihedral, camber, and twist of 
the aircraft wing. Thickness strips are used to model the 
thickness of the aircraft wing and pylons. [Ref. 10 :p. 807] 
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The capability exists to install multiple sets of wings, 
fins, or canards, and to use active control surfaces on the 
store by inputting the control laws into the program. Powered 
separations may be simulated by inputting the store thrust 
characteristics . 

The NEAR program actually consists of two separate 
programs: the source program and the trajectory program. Both 
are described in Section IV. B below. 

Alternate separation programs include USTORE and USSAERO 
codes. USSAERO was developed by F. A. Woodward, of NASA, as 
a lower-order panel method. USTORE was developed from USSAERO 
by G. J. van den Broek, of the National Institute for 
Aeronautics and Systems Technology, Pretoria, South Africa. 
[Ref. 11 : p . 309] 

B . THEORY 

The three principal tasks in the prediction of a store 
trajectory are: first, the determination of the nonuniform 
flow field in the neighborhood of the ejected store; second, 
the determination of the forces and moments on the store in 
this flow field; and third, the integration of the equations 
of motion to determine the store trajectory [Ref. 3]. 

1 . Source Program 

The source program is used to represent an 
axisymmetric body as a distribution of sources along the axis 
of the body. It provides point source— sink distributions to 
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represent the fuselage, rack, and store volumes. The program 
calculates and prints the source strengths and locations . 
These quantities are then used as input data to the trajectory 
program . 

The source program is used for the generation of an 
aircraft or pylon model which is then used for the trajectory 
program. These models are Mach number dependent and thus need 
to be generated for each test case with different Mach number. 
For ongoing store separation studies, a good database of 
aircraft models is required and should be available from 
appropriate Government facilities. The F/A-18 model was 
obtained from China Lake along with the program code. This 
model included pylon stations. 

2 . Trajectory Program 

The trajectory program uses the source-sink 
distributions from the source program, and additional 
information to first determine the vorticity distribution 
which represents the wing-pylon loading including interference 
of the fuselage, rack, and stores, and then to calculate the 
store trajectory. See Sections IV.C.3.d. (1) and (2) for 
descriptions of the vortex lattice method and the panel 
method . 

Once the trajectory program input has been generated 
for a single flight condition, it is relatively easy for the 
user to input different stores. The major effort is to obtain 
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the initial aircraft fuselage, wing, and pylon models for this 
input . 



3. Shape Modeling 
a . Aircraft 

The aircraft fuselage geometry was modeled by using 
a Fortran program called NGDELX developed by L. Gleason of the 
Naval Weapons Center. This program provides the coefficients 
to the NEAR polynomial representation of equivalent body. The 
shape is divided into appropriate segments. Radius values for 
the segment points and maximum radius values are entered. 
These coefficients are then entered into the source program to 
generate the source representation of the fuselage. Figure 9 
shows the coefficients for the F/A-18 fuselage. Following is 
the equation used to calculate the coefficients: 

r. q/L = C x + C 7 *(C 2 *(X/L) 2 + C 3 *(X/L)+C 4 ) + C s * (X/L) + C 6 * (X/L) 2 



Jb. Pylons 

The pylons used are modeled the same as in Section 
IV.C.l.a above. The geometry is used to obtain the NEAR 
coefficients . 



28 



Polynomial Rapraa tntation of Equivalent Body 
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Figure 9 F/A18 NEAR Polynomials 

c. Pod or Missile 

Again, the procedure in Section IV.C.l.a is used to 
represent the shape as a set of polynomial coefficients. The 
NGDELX program is again used to represent the missile shape as 
an "equivalent body of revolution" (EBR) . NGDELX is very easy 
to use, but, again, only represents axisymmetric shapes. 

d. Wing (Lifting Sv-fs.ces) 

The aircraft wings, missile fins, or any other 
lifting surfaces are modeled using a vortex lattice method to 
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represent the loaded wing. The thickness is also modeled 
using the lattice method. 

(1) Vortex Lattice Methods . There are several 
variations of the vortex lattice method that are presently 
available and have proven to be very practical and versatile 
theoretical tools for the aerodynamic analysis and design of 
planar and non-planar configurations. [Ref. 12 p. 27] The 
vortex lattice method represents the wing as a planar surface 
on which a grid of horseshoe vortices is superimposed. The 
velocities induced by each horseshoe vortex at a specified 
control point are calculated using the law of Biot-Savart . A 
summation is performed for all control points on the wing to 
produce a set of linear algebraic equations for the horseshoe 
vortex strengths that satisfy the boundary condition of no 
flow through the wing. The vortex strengths are related to 
the circulation and the pressure differential between the 
upper and lower wing surfaces. The pressure differentials are 
integrated to yield the total forces and moments. [Ref. 13 
p.261] For a rigorous introduction to the vortex lattice 
method, see Reference 13. 

(2) Panel Methods. The configuration is 
modeled by a large number of elementary quadrilateral panels 
lying either on the actual aircraft surface, or on some mean 
surface, or on a combination thereof. To each elementary 
panel, there is attached one or more types of singularity 
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distributions, such as sources , vortices, and doublets. 
These singularities are determined by specifying some 
functional variation across the panel (e.g., constant, linear, 
quadratic, etc.), whose actual value is set by corresponding 
strength parameters . These strength parameters are determined 
by solving the appropriate boundary condition equations . Once 
the singularity strengths have been determined, the velocity 
field and the pressure field can be computed. [Ref. 13:p. 258- 
259] 

C . PROCEDURE 

1 . Source Program 
a . Input 

The program input consists of the polynomial 
representation of the equivalent body of revolution (EBR) 
values obtained in Section IV.B.3.a. The surface of the EBR 
is approximated by these polynomials, which represent the EBR 
x, r distribution. The source program is then run with a user- 
specified finite distribution. The surface obtained from this 
source distribution can then be compared with the polynomial 
surface. Additional runs may be required to closely match the 
surfaces. The aircraft fuselage, pylon (s), and the adjacent 
store EBR values are input along with two variables, NRAT and 
PERCR. NRAT is the number of segments which the body will be 
divided for the specification of the source distribution. 
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PERCR is the source spacing for each NRAT segment and is input 
as a fraction of the local body radius of the segment. 

Adjacent stores are included in the aircraft source 
input. In this way it becomes part of the aircraft 
configuration. The separated store is not input to the source 
program since the trajectory program calls the source program 
during its execution to model the store. this is due to the 
fact the store changes position during the program run and the 
sources/sinks will change values. Figure 10 lists a sample 
format for the source program input. 
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Figure 10 Source Program Input 
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b. Output 

The source program output is used as input to the 
trajectory program. However, it is not directly read into the 
trajectory input and must be entered via keyboard. The form 
of the output is shown as a partial output in Figure 11. 
Source locations and strengths are listed. 
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Figure 11 Source Program Output 
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2 . Trajectory Program 

The trajectory program uses the output of the source 
program, and the input of the wing configuration and pylons, 
along with store information to calculate the trajectory of 
the separated store. It is not the intent of this report to 
fully explain the aircraft modeling details of the program. 
The assumption is made that the appropriate aircraft model is 
available for the correct flight condition and configuration. 
The applicable store inputs to the program are discussed and 
the resulting trajectory examined. Nielson contains a 
thorough description of the theory and input of the aircraft 
model in Reference 3 . 

a . Input 

The example input to the trajectory program was 
1064 lines in size. The store data entry begins at 
approximately 980 lines, depending upon aircraft input. The 
store input will vary according to the physical configuration 
of the store. For the example separation. Input Items 35-43 
and 48-51 were used to input store physical data. The ejector 
information was input in Items 44A-47. This is where the 
fifth-order polynomial representation of the ejector force was 
used. Data concerning initial and final times, and 
integration step sizes were entered in Input Item 72. The 
correct input form of the above data is contained in the NEAR 
code users manual [Ref. 14] and is not presented here. 
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b . Output 

The output of the trajectory program is dependent 
upon the input conditions, i.e., the number of stores, number 
of store segments, and number of integration steps, etc. The 
example case output is 4,524 lines. The first 1,840 lines of 
output are aircraft-specific output, such as the source and 
vorticity information. 

The last 2,650 lines of code contain the trajectory 
data. The ejector force data is repeated as is the store 
shape data, reference dimensions and inertial characteristics. 

Proceeding the ejector data, a data block is 
presented for each time integration step. This data block 
contains the parameters which describe the store displacement, 
velocity, and acceleration during the separation. A partial 
sample output of this section of the output file is included 
in Figure 12. Plotting routines can be developed to 
facilitate data analysis of the simulation but have not been 
developed for this study. The following are short 
descriptions of some of these parameters. A comparison with 
the results of the method of Section III follows in Section 
IV. D. 

The force and moment coefficients are listed by 
their individual contributions. The effects of buoyancy, 
slender body theory, crossflow, and empennage are totaled to 
give an effective C H , C Y , C^, and C ln . 
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Figure 12 Trajectory Program Fartial Output 
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Figure 13 Trajectory Output (Cont) 
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The ejector force and moment components are listed. 
The vertical force history matches the force profile from 
Section II very closely. The store's nose, inertial reference 
center, and the base position in the fuselage coordinate 
system is listed. The separation distance from the initial 
position is also listed as output. Plots of these parameters 
will indicate the miss distance of the shape from the wing. 
The store translational and rotational velocities are output 
in both the store' s reference frame and the fuselage reference 
frame . 
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V. 



RESULTS AND COMPARISONS 



A EJECTOR FORCES AND MOMENTS 

1 . Vertical Velocity 

The predicted end-of-stroke vertical velocity from the 
hand-calculation was 29.9 fps, comparing with the NEAR code 
results of 30.84. This is a 3% difference. The predicted 
value from the ejector graph of Reference 4 yields 24.5 fps. 
The ejector graph is representative of the average store used 
in the Naval inventory, which more than likely would have 
control surfaces for stability reasons. These surfaces would 
also provide aerodynamic damping which could be the source of 
the discrepancy. 

2 . Pitch Rate Q 

The predicted pitch rate, q, of 1.12 radians/second 
matches well with the NEAR end-of-stroke q of 1.08 
radians/second. For most applications, this correlation of 4% 
should have sufficient accuracy. 

B. TRAJECTORY COMPARISON 

The main consideration in determining the accuracy of 
these trajectory methods is the safe separation (or jettison) 
of the pod. Vertical distance and pitch attitude are 
indications of the pod separation while velocity and angular 
velocity are indications of the store loads. Therefore, 
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criteria for judging the correlation between the linear method 
and the NEAR code must be based on the practical standpoint of 
a reasonable assumption of safe separation. The individual 
parameters are addressed below. For determining the valid 
time range of the aerodynamic prediction method, the NEAR code 
was used as reference. 

1 . Vertical Separation Z 

Vertical distance z in the linear method is derived 
from the vertical velocity term. Therefore the error 
(deviation from NEAR results) will be slightly less due to the 
integration effects. However, given the initial vertical 
velocity, the NEAR code predicts a five foot separation. The 
error limit must be set as an acceptable variation of this 
distance. Based upon practical knowledge of store separation, 
a 20% error limit was set. Figure 14 compares the linear 
aerodynamic method results comparison with the NEAR results in 
the amplified time region of the first 0.2 seconds. 

2 . Pitch Angle 

Since the pitch angle is used solely for determining 
the safe separation and is not used for more accurate 
applications, like rocket engine ignition, the effect of an 
error in pitch will vary as the sine of the angle deviation. 
This could result in an error, due to rotation about the 
center of gravity, of an order of 7.2 feet times the sine of 
the error angle. The pitch angle error must be fairly small 
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initially, but can grow as time increases due to increased 
separation distance and its decreasing effect. Therefore, a 
conservative error limit is set for 0.1 sec. At 0.1 seconds, 
the NEAR pod has rotated 4.5 degrees. A 30% error would give 
a vertical error of 0.2 feet, which is approximately 10% of 
the actual separation. Figure 15 shows a comparison of the 
linear method with the NEAR results. Figure 16 provides the 
error between the methods as a function of time. 
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Figure 15 Pitch Angle 
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Figure 16 Parameter Error 
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TIME (SECS) 



3 . Pitch Rate 



The required accuracy for pitch rate will be 
determined by its application. For instance, navigational 
considerations will require quite accurate predictions of the 
angular loading of its components . However, for this case, 
there is not a requirement of this kind. The divergent nature 
of the pitch rate from the linear aerodynamic model does not 
present a problem for this case, beyond the fact that the 
pitch angle follows this trend. Therefore, an error limit is 
not placed upon the pitch rate. 

4 . Vertical Velocity 

Vertical velocity is divergent and represents the lack 
of aerodynamic damping in the linear model. Therefore, it 
does not provide meaningful information beyond the phase 
relationship with the vertical separation distance. For this 
reason, an error limit is not placed upon the velocity. 

C . AUGMENTED AEROPRED ICT ION 

Due to the limited valid time limits on the aerodynamic 
method, an alternative approach was investigated. Obviously, 
the non-linear operating region is reached very quickly. The 
two dominant terms of the A-matrix are the damping terms, C,*, 
and 0^. The unstable pod is shown from the NEAR code to 
quickly reach a somewhat steady-state value of pitch rate as 
it tumbles. An effort was made to calculate a new (and valid) 
value of C,^. This approach modeled the tumbling store as a 
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cylinder of constant radius. Using a value of C d =1.0, and the 
end-of-stroke pitch rate of 1.12 radians/second, new values of 
damping coefficients were calculated. 

1 . Ci*,- Coefficient 

The new pitch damping coefficient of -.1054 was used 
versus the MISDATCOM value of -.04. This is an obvious 
improvement and also realistic since the pod will definitely 
experience aerodynamic damping, due to pitch rate, that would 
not have been predicted by MISDATCOM. 

2. C**- Coefficient 

The value of the aerodynamic damping effects of change 
in angle of attack, C^, was combined with C^ in the MISDATCOM 
output. Therefore, the valid value of C^ is unclear. For 
this model of a rotating pod, it was assumed that C^ would 
have minimal effect on the damping and a value of 0.0 was 
assigned. 

3 . Trajectory Prediction 

The predicted parameters using the above coefficients 
are contained in Figures 17, 18 and 19. As shown in the 
figures, there is considerable improvement in the correlation 
with the NEAR trajectories for pitch angle and pitch rate. 
The more realistic value for C^ obviously improves the 
results. The time range of validity for these two parameters 
is now estimated to be 1.0 seconds. However, the most 
important parameter, vertical separation distance, is still 
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outside of reasonable accuracy limits. Therefore, the 
augmented plant matrix method is still unusable for complete 
trajectory prediction. Figure 20 shows the parameter error 
for the augmented system. 

Although there is substantial correspondence in the 
partial comparison between the NEAR results and these results, 
there are not enough data points to substantiate this method. 
Further work using different shapes need to be done. 
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Figure 18 Pitch Angle (Augmented) 
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Figure 19 Pitch Rate (Augmented) 
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Figure 20 Parameter Error (Augmented) 
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VI . CONCLUSIONS 



In the preliminary design application, MISDATCOM provides 
an easy-to-use method to predict moderately accurate store 
aerodynamic characteristics. These aerodynamic coefficients 
can be used to model the store dynamics for design work. 
However, their use in the linearized equations of motion for 
an externally ejected store trajectory prediction is highly 
limited . 

Comparisons with the Nielson Engineering and Research 
(NEAR) code show that trajectory values are outside of 
reasonable accuracy limits for use in store separations. The 
store separation vertical distance and pitch angle values 
diverged within 0.1 seconds. Therefore, the use of the 
linearized equations of motion cannot provide useful 
trajectory information. 

However, a modification to these coefficients to more 
accurately reflect the existing aerodynamic damping does 
provide correlations with the NEAR code of within 30% for 
pitch angle and 10% for pitch rate. Use of these "augmented" 
linear equations of motion for attitude prediction shows 
sufficient accuracy for many applications. 

The NEAR code is widely used for ejected store trajectory 
prediction in the aviation field. Its use for safe separation 
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investigation should be considered as standard for obtaining 
accurate results. The cost of the Nielson code in both 
engineering time and computational time increases with the 
level of information produced, but the need for safe 
separation accuracy, aerodynamic load distributions, and 
vortex-induced effects make the extra effort justified. 

For limited applications requiring less accuracy, the 
augmented linear equations of motion can provide some 
information for the initial trajectory movement. A potential 
application would be the investigation of the effects of small 
weight or inertial modifications to existing stores. 



51 



APPENDIX A 



A. POD CHARACTERISTICS 

The pod in question was an ALE-2 9B pod modified to 
accommodate the SUU-53 chaff dispenser and associated 
equipment. The ALE-29 pod was flight approved by the Aircraft 
Configuration Control Board (ACCB) document #75-004. It was 
originally an ALE-2 pod. No significant aerodynamic changes 
have been made to the pod. 

The pod physical characteristics are listed in Table 1 . 

TABLE 1 POD CHARACTERISTICS 



Length 


153.8 inches 


Radius -minimum 


9.75 inches 


-maximum 


11.20 inches 


CG location -empty 


62.4 inches from nose 


-full 


65 . 4 inches 


Weight -empty 


326 lbs 


-full 


371 lbs 


Moment of Inertias 




-Ixx 


300 lb-ft 2 


-iyy 


1675 lb-ft 2 


-Izz 


1675 lb-ft 2 
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APPENDIX B 



A. FORCE EQUATIONS 

IF = M* V = d/dt(M*V) = G 
Separating variables: 

2F x = G x lF y = G y lF z = G z 

Assuming lF x>y = 0, 

Integrating the remaining equation, 

|lF z dt = G 2 - Gj = (M*V z ) 2 - (M*V Z ) ! 

Assuming straight and level flight, V zl = 0 yields: 

JlF z dt = M* AV Z = M*V Z 

Integrating the ejector force polynomial yields a total 
impulse of 345 lb f -secs. Substituting (wt=371 lb), 

345 lb f _ secs = 371 lb m * V z 

V z = 1.1773 ( lb f -sec/ lbj * 32.174 ( lb m -ft / lb f -sec ? ) 

V z = 29.9 fps 
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B. MOMENT EQUATIONS 



ZM 0 = r x ZF = r x M*V = d/dt (r x M*V) = H 0 

where H 0 is the angular momentum about pt 0. Integrating the 
moments, 

/Zm o dt = H + r x H 

Separating the component equations yields, 



IM X = 


Ixx*w x - 


(Iyy - 


Izz) *W y *0) z 


II 

£ 


lyy*0)y ~ 


dzz " 


I xx ) * 0) z * 0) x 


II 

N 

£ 


Izz*W z - 


dxx - 


Iyy) *(O x *a)y 



Assuming the reference frame coincides with the principal 
axes, I xyfX2 , yz = 0, and assuming cl) X jZ = 0, yields the following: 

/lM y dt = I yy * My 

(Oy =1.12 rad/sec 
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APPENDIX C 



A. LONGITUDINAL EQUATIONS OF MOTION 

The equations of motion for the pod can be derived from 
Newton's Second Law of motion, which states that the summation 
of all external forces acting on a body must be equal to the 
time rate of change of the momentum of the body, and the 
summation of the external moments acting on a body must be 
equal to the time rate of change of the moment of momentum 
(angular momentum) . The time rates of change are all taken 
with respect to body coordinates space. [Ref. 8] 
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£ C^F x -m* ( U+W*Q-V*R) 

£ &F z =m* ( W+ V*P~U*Q) 

£ A M=Q*I y + P*R* (I x -I z ) + ( P 2 -R 2 ) *J XZ 



&F x =m* ( U+W*Q-V*R) 

£ A F z =m* (W+V*P-U*Q) 

£ &M=Q*I y 

By restricting the disturbances to small perturbations 
about the equilibrium condition, the product of the variations 
will be small in comparison with the variations and can be 
neglected, and the small angle assumptions can be made 
relative to the angles between the equilibrium and disturbed 
axes . 



£ AF x =m*ii 

&F z =m* (w-U*q) =m* (w-C/*0) 
£ AM=I y *<j=I y *6 



Expanding the applied forces and moments in terms of the 
changes in the aerodynamic and gravitation forces and moments 
using the total differential form yields: 
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£ dF x ~ (dF x /dU) *dU+ ( dF x /dW ) *dw+ ( dF x /dW ) *dW+ (dF x /d6) *dS- (dF x /dS, 
dF 2 = ( dFjdU ) *dU+ [dFjdW) *dw+ (d F z /dw) *dw+ (dF z /d 6) *d9~ (8 f 2 /60 
r AW= {dM/dU) *dU+ (dM/dW) +dW+ (dM/dW) *dW+ (dM/d&) *d& 



Non-dimensionalizing the terms in the above equations 
yield the aerodynamic coefficients. Following are the 
definitions of the longitudinal stability coefficients and 
derivatives : 

C IU Variation of drag and thrust with u. 

C xa Lift and drag variations along the x-axis . 

C v Gravity 

C I( i Downwash lag on drag. 

C xq Effect of pitch rate on drag. 

C IU Variation of normal force with u. 

C ia Slope of the normal force curve. 

C zii Downwash lag on lift of tail. 

C iq Effect of pitch rate on lift. 

C mu Effects of thrust, slipstream, and flexibility. 

C mct Static longitudinal stability. 

C mix Downwash lag on moment. 

C mq Damping in pitch. 

Arranging the linearized longitudinal equations into 
matrix form yields the following matrix equation. The contx*ol 
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matrix B is not shown here but is included in the Control-C 



program. Since there are no active control surfaces, the B- 
matrix is a zero matrix. 



Au 




x u 




0 


-g 




"Au 


Aw 




Zu 


z„ 


u 


0 




Aw 


Aq 


= 


M U +M„*Z U 


Mw+M»*Z w 


Mq+M^U 


0 




Aq 


E 

> 

<x> 




0 


0 


0 


0_ 




.A9, 
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B. LATERAL EQUATIONS OF MOTION 

In the same manner used to obtain the longitudinal 
equations, the lateral equations are derived. The same 
assumptions using perturbation theory apply. Following are 
the equations in their developing order. 



Y A F y =m* ( V+U*R-W*P) 

Y kL=P*I x ~R*J xz +Q*R* (I z -I y ) -P*Q*J XZ 

Y &N=R*I Z +P*J XZ +P*Q* (I y -I x ) + Q*R*J XZ 

Y A Fy-m* ( v+U*r + u*r) 

Y AL=£>*l x -r*J xz 
Y A N=r*I z -£*J xz 

£ AF y -w* ( v+CJ*r) 
y k.L=p* I x - r*J xz 



Y dF y ~ (dF y /d B) *dB + ( dF y /d Y) *cTP + (d F y /d<b) *dfl> + (d F y /84>) *d4>- ( dF y /- 



Apr 




Yp/u Y p /u — ( 1— Y r /u) g*cos0/u 




A (3 


Ap 




L p L p L r 0 




Ap 


Ar 


= 


Np N p N r 0 




Ar 


A<t> 




0 1 0 0 

L J 




<i 

i 



59 



APPENDIX D 



m-input( 'entet m ( 1 bs ) ! 9 ) /3 2 . 2 , 

Uo-input( 'enter Uo (fps)j'); 

S-input( 'enter S ( 1 n * * 2 ) t 9 ) /I 4 4 , 
b- i npu t ( ' en te r b (in)i'}/12; 

Ro- i npu t ( ' en t e r Ro ( b lugs/ft ** 3 ) i ’) f 
c r e f - 1 npu t ( ' en te r c r e f ( i n ) t ’ )/l 2 , 

Ix-lnput( ' enter Ixx ( f t - 1 b ) t ' ) /3 2 . 2 / 

Iy-input( 'enter IYY ( f t- 1 b ) : ' ) /3 2 . 2 ; 

Iz-lnput( 'enter Izz ( f t- lb ) s ' )/3 2 . 2 , 

9-32 . 1 74 | 

Cn-lnputj 'enter Cn (Normal Force coe f f i c 1 e n t ) ! ' ) ; 

Cm-input( 'enter Cm (Pitching Moment Coe f f 1 cl en t ) s ’ ) f 
Ca-inputj 'enter Ca (Axial Force Coe £ f i c i e n t ) i ' ) , 

Cy- i npu t ( ' en te r Cy (Side Force Coe f £ 1 ci en t ) t ' ) ; 

CLN-input( 'enter CLN (Yawing Moment Coeff):')> 

CLL-lnput( 'enter CLL (Rolling Moment Coeff)i')j 
CNA-input( 'enter CNA (Normal Force Deriv wrt AOA ) : ' ) ; 

CMA-input( 'enter CM A (Pitching Moment Deriv wrt AOA ) t ' ) , 

CYB-input( 'enter CYB (Side Force Deriv wrt Beta):'), 

CLNB- i npu t ( 9 e n t e r CLNB (yawing Moment deriv wrt Beta),'), 
CLLB-lnput( 'enter CLLB (Rolling Moment deriv wrt Beta)i'), 

CL- 1 npu t ( ' ente r Cl (lift coeff)i')j 
CD- input (' ente r Cd (Drag coeff)i'), 

CNQ- i npu t (' ente r CNQ (Normal Force Deriv wrt Q ) j r ) | 

CMQ- 1 npu t (' ente r CMQ (Pitching Moment Deriv wrt Q)i'), 

CNAOADOT- i npu t (' ente r CN AOAdot (Normal Force Deriv wrt AOAdot)! '), 
CMAOADOT-input (' enter CM AOAdot (Pitching Mom Deriv wrt AOAdot)!'), 
q- . 5 * ro *uo *uo , 

Xu--2*Cd*q*s/(m*uo) , 

Xw-Cl *q*8/(m*uo ) , 

Zu--2*Cl *q*8/(m*uo) , 

Zw--Cd*q*8/(m*uo) , 

Zwd-Cnaoadot*cref*s*q/(2*uo*uo*m), 

Mw-Cma *q* s *cre f/( Iy *uo ) , 

Mu-0 | 

Ma-uo *Hw| 

Mwd-Cmaoadot*cref*q*i*cr«f/( 2*uo*uo*Iy) , 

Mad-uo *Mwd | 

Mq-Cmq*cref *q*8*cref/( 2*uo* Iy ) , 

A-(Xu Xw 0 -q;Zu Zw uo 0 , ( Hut Mwd * Zu ) /( 1 -Zwd ) ( Mw+Mwd * Zw ) / ( 1 - Zwd ) . . . 

( Mq+Mwd*uo )/( 1-Zwd ) 0|0 0 1 OJ; 

B - ( 0 0 0 0,0 0 0 0 J ' I C- 1 0 1 0 0 | | D- ( 0 Oj, 

IC-( 0 0 0 0 ) ' I 

IC( 2 )-lnput( ' enter Initial vertical velocity (fps)i'), 

IC( 3 ) - 1 npu t ( ' ente r initial pitching rate (rad/s)i'), 

t t t- 1 npu t (' ente r end time (sec)')) 

t-0i .OOlittt, 

uuu-( ttt/. 001 ) 4 1 , 

u-0 # ones( 2,uuu) , 

8 i mu ( ' IC' , IC) , 
y-simu(a,b,c,d,u,t) , 
erase , 
plot( t,y ) I 



Figure 22 Control— C Program For EOM Simulation 
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Figure 23 MISDATCOM INPUT/OUTPUT 
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Figure 24 
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NOTE - A INDICATES SLOPE DISCONTINUOUS POINTS 
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Figure 26 
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Figure 27 
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